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Abstract 

In  this  paper,  we  report  on  a  new  stability  analysis  for 
hybrid  legged  locomotion  systems  based  on  factoriza¬ 
tion  of  return  maps.  We  apply  this  analysis  to  a  fam¬ 
ily  of  models  of  the  Spring  Loaded  Inverted  Pendulum 
(SLIP)  with  different  leg  recirculation  strategies.  We 
obtain  a  necessary  condition  for  the  asymptotic  sta¬ 
bility  of  those  models,  which  is  formulated  as  an  ex¬ 
act  algebraic  expression  despite  the  non-integrability 
of  the  SLIP  dynamics.  We  outline  the  application  of 
this  analysis  to  other  models  of  legged  locomotion  and 
its  importance  for  the  stability  of  legged  robots  and  an¬ 
imals. 


1  Introduction 

This  paper  introduces  a  new  formalism  for  studying 
the  stability  of  legged  locomotion  gaits  and  other  pe¬ 
riodic  dynamically  dexterous  robotic  tasks.  We  are 
most  immediately  motivated  by  the  need  to  explain 
and  control  the  remarkable  performance  of  RHex,  an 
autonomous  hexapedal  robot  runner  with  unparalleled 
mobility  [1].  Powered  by  only  six  actuators,  located 
at  the  “hips”  to  drive  each  of  its  six  passively  com¬ 
pliant  legs,  RHex’s  locomotion  is  excited  by  a  single 
periodic  “clock”  signal  split  into  phase  and  anti-phase 
copies  for  coordinating  its  alternating  tripod  gait.  A 
simple  PD  controller  at  each  hip  motor  in  a  given  tri¬ 
pod  forces  its  leg  to  track  the  alternately  fast  and  slow 
clock  reference  signal  corresponding  to  putative  stance 
and  swing  phases.  Although  RHex  is  currently  being 
endowed  with  more  sophisticated  proprioceptive  reflex 
loops  to  increase  its  responsiveness  to  varied  terrain, 
we  are  strongly  motivated  to  develop  an  analytical  un¬ 
derstanding  of  the  relationship  between  clock  signal 
and  steady  state  gait  properties  in  the  “simple”  open 
loop  case. 

A  complete  account  of  this  relationship  in  even  the 
simple  case  would  entail  insight  into  the  steady  state 


properties  of  an  under-actuated  high  degree  of  freedom 
hybrid  mechanical  system  whose  Lagrangian  dynam¬ 
ics  switches  between  a  set  of  2®  possible  holonomically 
constrained  models  depending  upon  which  toes  are  in 
contact  with  the  ground.  Fortunately,  a  growing  body 
simulation  study  and  empirical  evidence  [2]  suggests 
that  RHex,  when  properly  tuned,  exhibits  sagittal 
plane  stance  behavior  well  approximated  by  the  two 
degree  of  freedom  SLIP.  Thus,  in  the  short  term,  we 
seek  to  understand  how  adjustments  to  a  coordinat¬ 
ing  clock  signal  will  determine  the  steady  state  perfor¬ 
mance  of  the  bipedal  SLIP  underlying  the  alternating 
tripod  gait.  In  this  paper,  we  develop  the  mathemat¬ 
ical  foundations  of  a  new  formalism  for  distinguishing 
volume-preserving  from  non- volume-preserving  hybrid 
Lagrangian  systems.  We  apply  this  formalism  to  the 
hybrid  SLIP  model,  the  results  of  which  suggest  new 
insight  into  the  relationship  between  clock  excitation 
and  steady  state  gait. 


1.1  Leg  Swing  Policies  and  Self  Stabil¬ 
ity  in  the  SLIP  Template 

The  SLIP  model  provides  a  ubiquitous  description  of 
biological  runners  in  the  sagittal  plane  [3]  and,  as  men¬ 
tioned  above,  a  broadly  useful  prescription  for  legged 
robot  runners  such  as  RHex  [4,  1,  2]  as  well.  The 
closely  related  three  degree  of  freedom  Lateral  Leg 
Spring  (LLS),  has  been  recently  identified  as  a  candi¬ 
date  template  for  cockroach  running  in  the  horizon¬ 
tal  plane  [5,  6]  and  seems  likely  to  be  relevant  for 
RHex  as  well  [1].  For  present  purposes,  the  most  im¬ 
portant  insight  from  these  models  has  been  to  pro¬ 
vide  a  mathematical  explanation  for  their  unexpected 
“self-stability”  properties  (asymptotically  stable  equi¬ 
librium  gaits  in  the  absence  of  any  sensor  based  feed¬ 
back  inputs),  thus  affording  new  hypotheses  regarding 
biological  control  strategies  [7]  and  a  formal  founda¬ 
tion  for  elucidating  the  stable  open  loop  performance 
of  RHex. 
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The  originally  discovered  self-stability  of  SLIP  [8,  9]  is 
based  on  a  simple  leg  swing  policy,  the  leg  touchdown 
angle  at  the  end  of  aerial  phases  being  kept  constant 
from  stride  to  stride.  A  similar  strategy  was  earlier 
shown  to  yield  stable  gaits  in  the  LLS  model  [6].  Re¬ 
cently,  a  different  time-dependent  leg  retraction  pol¬ 
icy  has  been  shown  to  inherit  the  stability  properties 
of  the  fixed  touchdown  angle  policy  while  increasing 
the  robustness  of  the  SLIP  system  [10].  On  the  other 
hand,  numerical  simulations  of  a  recirculation  policy 
where  the  SLIP’s  leg  starts  recirculating  after  leg  liftoff 
at  a  constant  angular  velocity  until  leg  touchdown  sug¬ 
gest  not  asymptotic  but  neutral  stability  [11].  Hence 
the  leg  swing  policy  seems  to  play  a  central  role  in  the 
stability  of  those  low-dimensional  models,  as  is  also 
suggested  by  high  dimensional  systems  such  as  RHex 
[1]  or,  as  established  in  an  independently  conceived 
model  inspired  by  animal  locomotion  strategies,  for  a 
quadrupedal  trotter  [12]. 

Stability  of  such  systems  is  observed  through  the  full- 
stride  return  map  -  the  function  relating  body  state  at 
one  stride  to  body  state  at  next  -  which  summarizes  all 
properties  relevant  to  the  goal  of  translating  the  body 
center  of  mass.  Unfortunately,  even  the  simplest  2dof 
SLIP  system  is  non-integrable  [13],  which  precludes 
closed  form  solutions  that  might  illuminate  stability 
mechanisms.  In  this  paper  we  will  show  how  the  sta¬ 
bility  of  hybrid  systems  possessing  certain  symmetries 
can  be  analyzed  in  terms  of  their  non-hybrid  compo¬ 
nents,  e.g.  flight  phase  and  stance  phase  in  the  case  of 
SLIP,  thus  decomposing  the  hybrid  return  map  into 
“partial”  return  maps  that  might  be  analyzed  more 
easily. 

Before  doing  so,  we  first  introduce  certain  terminology 
and  notation  used  throughout  the  paper  by  way  of 
reviewing  Liouville’s  theorem  (see  e.g.  [14]),  pointing 
out  that  the  classical  result  should  not  be  expected  to 
apply  in  the  present  case. 

1.2  Liouville’s  theorem  and  stability 

Liouville’s  theorem  states  that  volume  in  phase  space 
of  a  holonomically  constrained  conservative  dynami¬ 
cal  system  described  by  a  single  Hamiltonian  flow  is 
preserved,  i.e.  a  set  of  initial  conditions  at  t  =  to  in 
phase  space  will  be  mapped  to  a  set  with  identical 
symplectic  volume  for  any  t  >  to- 
In  the  case  of  non-holonomically  constrained  conser¬ 
vative  Hamiltonian  systems,  the  number  of  indepen¬ 
dent  conjugate  momenta  is  smaller  than  the  number 
of  independent  configuration  space  variables  and  Liou¬ 
ville’s  theorem  cannot  be  applied  (see  e.g.  the  asymp¬ 
totic  stability  of  the  Chaplygin  sleigh  in  [15]  and  ref¬ 
erences  therein). 


In  the  case  of  piecewise-defined  holonomically  con¬ 
strained  conservative  Hamiltonian  systems  with  dif¬ 
ferent  flows  with  a  in  an  indexing  set  /  the  tran¬ 
sition  to  a  new  Hamiltonian  flow  is  triggered  by  so- 
called  threshold  functions  (for  a  general  definition  of 
hybrid  systems  see  [16]).  In  almost  all  settings  within 
robotics,  these  threshold  functions  between  different 
flows  /*  and  /|  depend  upon  state,  and  often  have  no 
explicit  time  dependence  at  all.  Examples  include  a 
discrete  version  of  the  Chaplygin  sleigh  [15,  17]  and 
low-dimensional  models  of  legged  locomotion  in  the 
horizontal  and  sagittal  plane  [6,  8,  9],  which  all  ex¬ 
hibit  partial  asymptotic  stability  for  certain  param¬ 
eter  settings.  Here,  the  (local)  asymptotic  stability 
of  those  hybrid  system  at  a  fixed  point  means  that 
the  eigenvalues  of  the  components  of  their  linearized 
return  maps  defined  by  Poincare  sections  lie  within 
the  unit  circle  (excluding  those  corresponding  to  con¬ 
served  quantities).  Liouville’s  theorem  is  not  directly 
applicable  due  to  the  hybrid  nature  of  such  systems, 
despite  the  fact  that  the  flows  are  Hamiltonian. 

In  those  cases,  the  return  map  TZ  defined  by  a  Poincare 
section  is  composed  of  several  “factor”  maps  that 
relate  the  state  variables  directly  after  one  transition 
to  those  directly  before  the  next.  Additional  transition 
mappings  are  used  at  transitions  from  flows  to 
/^.  They  arise  in  practice  from  the  fact  that  different 
flows  are  most  easily  handled  analytically  in  different 
coordinate  systems.  Thus  TZ  =  o  r“  o  •  •  •  o  o  r^. 

1.3  Contribution  of  this  paper 

In  this  paper,  we  focus  on  the  role  of  volume  preserva¬ 
tion  in  flows  and  transition  maps  of  models  of  legged 
locomotion  as  an  indicator  of  local  stability.  Although 
we  cannot  invoke  Liouville’s  theorem  for  reasons  just 
reviewed,  we  nevertheless  deduce  the  necessity  for  vol¬ 
ume  preservation  at  a  given  fixed  point  in  these  flows, 
manifesting  itself  as  the  condition  of  unity  determi¬ 
nant  in  the  associated  linearized  return  map.  In  gen¬ 
eral,  in  order  to  check  whether  the  return  map  of  a  hy¬ 
brid  system  TZ  is  volume-preserving  at  a  fixed  point, 
the  map  must  be  computed  explicitly  [17,  6,  8,  9]. 
However,  we  will  show  that  if  all  vector  fields  fa,  a  G  I 
possess  time  reversing  symmetries  and  if  the  peri¬ 
odic  orbit  giving  rise  to  a  fixed  point  of  TZ  is  composed 
of  pseudosymmetric  orbits  (to  be  defined  below)  on 
each  flow  domain  Va,  volume-preservation  of  the  whole 
system  can  be  determined  by  volume-preservation  on 
individual  flow  domains.  Moreover,  if  the  transition 
function,  h^,  from  a  flow  domain,  Va,  to  a  next  flow 
domain,  Vf},  enjoys  a  certain  symmetry  related  to  the 
vector  field’s  time  reversing  symmetry  Sa,  then  vol¬ 
ume  preservation  on  Va  can  be  determined  in  certain 
cases  without  an  explicit  expression  for  the  flow  map 
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on  that  domain. 

This  paper  introduces  the  rigorous  foundations  for  the 
formalism  just  outlined,  and  provides  an  example  of 
its  value  for  robotics  by  application  to  the  LLS  [6] 
and  SLIP  [8,  9]  models.  As  explained  above,  the  ques¬ 
tion  of  whether  a  hybrid  system  is  volume  preserving 
or  not  has  immediate  consequences  for  gait  stability. 
Since  our  conditions  for  volume  preservation  of  the  full 
stride  return  map  may  be  checked  with  respect  to  a 
simple  ’’factor  map”,  we  are  able  to  study  the  effects 
of  a  broad  range  of  leg  swing  policies  that  are  imple- 
mentable  in  the  analytically  tractable  ’’flight”  phase  of 
the  leg.  Specifically,  we  show  how  the  volume  preserv¬ 
ing  properties  of  SLIP  under  different  leg  recirculation 
strategies  can  be  determined  by  simple  differentiation. 
In  particular,  in  Sec. 3. 1.3  we  give  for  the  first  time 
necessary  conditions  for  asymptotic  and  neutral  sta¬ 
bility  of  a  RHex-like  leg  recirculation  scheme  applied 
to  the  SLIP  model,  without  explicit  computation  of 
the  stance  phase  flow  map. 

2  Factoring  return  maps  with 
time  reversing  symmetries 

Assume  a  hybrid  mechanical  system  whose  time  evo¬ 
lution  is  described  by  holonomically  constrained  au¬ 
tonomous  conservative  vector  fields  fa,  a  G  I  with 
configuration  space  variables  qa-  Xa  =  fa{xa)  with 
X  =  {qa  qa)^  G  Va-  The  open  flow  domains  Va  are 
called  charts  [16].  Transitions  between  vector  fields  are 
governed  by  threshold  functions  which  can  depend 
on  the  initial  condition  Xao  =  Xa{t  =  Q)  &  Va,  time 
t,  and  the  current  state  ff{xao)-^  We  demand  that 
for  each  chart  there  is  only  one  threshold  function  h^. 
Transitions  to  the  vector  field  fp  are  uniquely  defined 
by  fao(^ao)  =  mint>o{t  :  h^{ff{xao),XaO,t)  =  0}.  In 
addition,  all  transition  mappings  T£  are  assumed  to 
be  volume  preserving.  The  flow  map  for  the  ath 
vectorfield  is  implicitly  defined  via  t^Q  by  :  Xao 

tP 

fa‘°{xao),  where  Xao  is  assumed  to  be  the  result  of  a 
preceding  chart  transition. 

We  will  define  a  partial  return  map  for  evolutions  on 
individual  charts  Va-  Under  the  assumption  that  a 
periodic  orbit  of  the  whole  dynamical  system  is  com¬ 
posed  of  pseudo-symmetric  orbits  (to  be  defined  be¬ 
low)  on  individual  charts,  volume  preservation  of  the 
full  return  map  TZ  can  be  determined  in  some  cases 
without  explicitly  computing  r^. 

Assume  that  the  vector  field  fa  admits  a  linear  in- 


^Note  that  this  definition  generalizes  [16],  where  only 
depends  on  /‘ (xao)- 


volutive  time  reversing  symmetry^  Sa  {Sa  ■  Xa  = 
fa{Sa  ■  Xa)  Xa  =  fa{xa))  with  =  id.  This  im¬ 
plies  SaO  ffoSa  =  fa*.  We  now  investigate  a  compo¬ 
sition  of  two  partial  return  maps,  Xai  =  Sa  o  rd(xao), 
followed  by  Xa2  =  SaOrd(xai).  If  the  functional  iden¬ 
tity  =  tfg  holds,  then 

SaOr^oSaOrd  =  SaO  fa°‘°  O  Sa  O  fa"*’  = 

-t’” 

=  id{l) 

We  now  need  a  sufficient  condition  for  t^i  =  This 
is  given  by 

(hi{fa"*’{Xal),Xal,t1o)  =  0 

h'i{ff°{Xao),XaO,tio)  =  o)  ^  =  tfg  (2) 

which  can  be  rewritten  as 

{hiifa"*’  O  O  fa"*’{Xao),  Sa  O  fj*’ {x  ao)  ,  O  = 

h^{Sa{Xaflj,SaO  fa°{Xao),tao)  =  o) 
^  (hi{fa"°{Xao),XaO,tio)  =  o)  (3) 

which  can  be  considered  an  invariance  of  the  thresh¬ 
old  function  under  the  partial  return  map  Sa  o  r^  : 
XaO  Sa  ord(xao).  This  essentially  checks  that  r^  is 
“compatible”  with  time  reversibility  of  fa- 
Given  that 

5a  o  rf  o  5a  o  rf  =  id  (4) 

then 

D{Sa  orloSaO  rf)  = 
DSa-Dri{SaOri)-DSa-Dri  =  1  (5) 

Next,  we  call  a  trajectory  on  a  chart  Va  pseudosym- 
metric  if  a  fixed  point  Xa  of  the  partial  return  map 
Sa  o  I’a  with  Sa  o  rd{xa)  =  Xa  exists.  Evaluation  of 
expression  (5)  at  such  a  fixed  point  then  allows  us  to 
determine  the  square  of  the  determinant  of  the  Jaco¬ 
bian  of  rf. 

DSa  ■  Drd{xa)  ■  DSa  ■  Drd{xa)  =  1 

det{DSa  ■  Drd{xa)  ■  DSa  ■  Drd{xa))  =  1 

det"^  {Dr  d{x  a))  =  1  . 

Here  det^(iJ5Q)  =  1.  Hence  if  a  periodic  orbit  de¬ 
scribed  by  a  fixed  point  of  the  full  return  map  TZ  is 
composed  of  pseudosymmetric  orbits  on  Va  and  con¬ 
ditions  (2)  hold  on  each  chart,  the  full  return  map  TZ 

^The  importance  of  time  reversing  symmetries  for  the 
asymptotic  behavior  of  dynamical  systems  has  been  emphasized 
in  [18]. 
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is  volume-preserving  at  this  fixed  point.  If  conditions 
(2)  do  not  hold  or  if  periodic  orbits  of  the  full  return 
map  are  not  composed  of  pseudosymmetric  orbits,  no 
conclusions  can  be  drawn  from  this  argument.  Note 
that  here  “volume”  is  defined  with  respect  to  the  state 
variables  chosen,  and  in  general  is  not  a  phase  space 
volume,  and  that  the  computation  (at  the  fixed  point) 
is  local. 


3  Application  to  hybrid  models 
of  legged  locomotion 

3.1  SLIP  with  pitching 

The  SLIP  model  consists  of  a  rigid  body  of  mass  m 
and  moment  of  inertia  I  with  a  massless  springy  leg 
attached  to  an  unactuated  hip  joint  which  is  a  dis¬ 
tance  d  away  from  the  center  of  mass  (for  details  see 
[9]).  A  full  stride  consists  of  a  stance  phase  with  the 
foothold  fixed  and  the  leg  under  compression,  and  a 
flight  phase  in  which  the  body  describes  a  ballistic 
trajectory.  Hence  there  are  two  vector  fields  /i  (for 
stance)  and  /2  (for  flight)  and  the  return  map  can  be 
written  as  TZ  =  T2  o  o  o  r\.  We  assume  that 
a  periodic  orbit  of  period  one  is  composed  of  pseu¬ 
dosymmetric  stance  and  flight  phases.^ 


3.1.1  Stance  phase  of  SLIP  with  pitching 

The  equations  of  motion  that  describe  the  stance 
phase  of  SLIP  with  mass  m  and  moment  of  inertia 
/  read  in  the  conventions  of  Fig.  I  (see  also  [9]): 

C  =  —  g cos  ijj  —  ((  +  dcos  (tjj  +  6)) 

mr] 

..  .  .  f) 

=  —2(^ij}-\-gs\TL'il)  +  d— - (sin  (^/^ -I- 0)) 

mr] 

0  _  sin  _l_ 

rjl 

with  spring  potential  V{r])  where  the  compressed 
spring  length  77  =  d'^  +  +  2dC  cos  {ip  -I-  9)  <  770, 

the  relaxed  spring  length  being  denoted  by  rjQ.  The 
distance  between  the  ‘hip’  pivot  and  the  mass  center 
is  d  and  all  joints,  including  the  toe-ground  pivot,  are 
assumed  frictionless  and  moment-free. 

The  linear  involutive  time  reversing  symmetry  of 

(6)  acting  on  xi  =  {(  tp  9  (  tp  9)^  is  given  by 

5i=diag(l,-l,-l,-l,l,l)  .  (7) 

®This  was  proven  in  [9]  for  SLIP  without  pitching  dynamics 
and  without  gravity  in  stance. 


Figure  1:  Parametrization  of  SLIP  with  pitching  dy¬ 
namics. 


Under  5i,  the  spring  length  77  transforms  as  ^1(77)  =  77. 
Transitions  to  flight  occur  when  the  spring  length  has 
reached  its  rest  length  770  which  is  also  the  initial  spring 
length  at  the  beginning  of  stance.  Hence  the  threshold 
function  can  be  written  as  =  rf{t)  — 

tJq.  Then  the  partial  return  map  Si  o  rf  is  volume¬ 
preserving  at  a  fixed  point,  because  (3)  holds: 

hi{xi{tjo),xio,tjQ)  =r]^{tlo)-r]^  =  0 

^  hl{Si- Xio,Si- Xi{tl^),tlQ)  =r]l-r]'^{tlQ)  =  0 

This  result  is  independent  of  the  specific  form  of  the 
spring  potential  V{r]).  For  d  =  0,  the  two-dimensional 
SLIP  system  describing  a  point  mass  on  a  massless 
spring  is  recovered. 

3.1.2  Flight  phase  of  SLIP  with  pitching 

The  equation  of  motion  of  the  center  of  mass  that 
describe  the  flight  phase  of  SLIP  read 

y  =  0 

z  =  -g 

9  =  0  (8) 

where  the  z-axis  points  vertically  upwards,  the  77-axis 
points  in  a  horizontal  direction,  and  9  denotes  the 
SLIP  body’s  pitching  angle  with  respect  to  the  hori¬ 
zontal.  The  linear  involutive  time  reversing  symmetry 
S2  of  (8)  acting  on  X2  =  {y  z  9  y  z  9)^  is  given  by 

52  =  diag(±l,l,-l,Tl,-l,l)  .  (9) 

Here,  the  sign  ambiguity  in  9  was  resolved  by  match¬ 
ing  9’s  transformation  to  that  of  the  stance  phase  (7). 
The  transformation  law  of  y  under  ^2  is  not  needed  in 
subsequent  calculations  and  is  left  unresolved. 

The  simplicity  of  the  equations  of  motion  (8)  allows  us 
to  explicitly  compute  the  determinant  of  the  Jacobian 
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in  terms  of  partial  derivatives  of  4)(t) 


of  the  partial  return  map  S2  o  at  a  fixed  point  for 
a  given  leg  recirculation  scheme.  Therefore  the  appli¬ 
cation  of  the  formalism  of  Sec.  2,  which  only  provides 
a  sufficient  condition  for  volume-preservation,  but  not 
for  non-preservation,  is  relegated  to  the  appendix. 
The  equations  of  motion  for  the  2:  and  6  coordinates 
can  be  explicitly  solved  and  read  in  dimensionless  vari¬ 
ables: 

.  -  P 

Z{t)  =  ZQ+ZQt-  — 

z{i)  =  zo-i 

9{i)  =  Oq  +  Ogi 

ki)  =  k  (10) 


A  1  num 

det  {Drl{i2))  =  1  +  (13) 

with 

=  v^l-(5o  +  dcos(0o))"- 
(^~9k^4'{iTD)  +  ZQdzo(f){tTD)  +  ^odg^4’(iTDkj 
+zo  —  dsin{do)9o 

^iden  ^  -1^1  -  [So  +  d cos{0oWdf (pitro) 

-Zo  -I-  d  sin (00 )  6*0 


with  i  =  z  =  z  =  00  =  6*0,  and 

0  =  We  now  want  to  explore  different  strate¬ 

gies  to  position  the  leg  during  flight.  Since  the  leg  is 
assumed  massless,  any  leg  angle  trajectory  <p{i)  where 
(j)  is  defined  in  Fig.  1  can  be  commanded. 

The  threshold  function  for  a  recirculating  leg  reads 
in  dimensionless  variables 

hl{x2{i),X2o,i)  =  z{i)  +  dcos{9{i))  -  cos(^(t))  (11) 

where  d  =  ^.  Setting  eq.  (11)  to  zero  determines  the 
time  from  leg  liftoff  {Ilo  =  0)  to  leg  touchdown  Itd  ■= 
^20 1  for  which  in  general  a  closed  form  solution  does  not 
exist.  Then  the  flow  map  takes  the  dimensionless 

state  vector  X2  =  {z  9  z  9p  {y  and  y  are  omitted)  from 
its  value  at  leg  liftoff  to  that  at  touchdown:  r2(i2o)  = 
X2{tTD)-  A  fixed  point  of  a  pseudosymmetric  flight 
trajectory  satisfies  X2  =  S20  r2{x2)- 
The  determinant  of  the  Jacobian  of  can  easily  be 
computed  from  the  flight  trajectory  expressions  (10), 
bearing  in  mind  that  the  flight  time  Itd  also  depends 
on  the  initial  conditions: 

det(Ilr2)  =  1  —  d-J-TD  +  zodpTO  +  9odgJ,TD  (12) 

In  this  expression,  the  leading  term  1  is  a  consequence 
of  Liouville’s  theorem,  because  z  and  z  and  0  and  0  are 
canonically  conjugate  up  to  a  trivial  rescaling,  whereas 
the  remaining  terms  make  the  non-applicability  of  Li¬ 
ouville’s  theorem  to  this  hybrid  system  with  a  state- 
dependent  threshold  function  (11)  explicit.  Albeit  Itd 
cannot  be  computed  explicitly  in  general,  tro  =  2zo 
at  a  fixed  point  of  S2  o  r\  and  cos(^(tTD))  =  zq  + 

(icos(0o),sin((/)(t7’£)))  =  —\Jl  —  (zq  +  (icos(0o))^  and 

9{tTD)  =  —9o.  Hence  using  implicit  differentiation  of 
(11)  the  determinant  at  this  fixed  point  can  be  written 


The  eigenvalues  of  the  partial  return  map  S2  o  are 
Ai  =  1  (vertical  energy),  A2  =  1  (rotational  energy), 
A3  =  —1,  and  A4  =  —  det(I?r2(i2)).  It  must  be  em¬ 
phasized,  however,  that  the  eigenvalues  of  this  partial 
return  map  are  not  equal  to  the  eigenvalues  of  the  to¬ 
tal  return  map  TZ  at  the  fixed  point.  In  particular,  the 
eigenvalues  of  TZ  can  be  complex,  as  in  Fig.  2,  below. 
Setting  eq.  (13)  to  1  yields  a  partial  differential  equa¬ 
tion  for  leg  recirculation  schemes  (j>{t)  that  are  vol¬ 
ume  preserving.  In  the  following,  for  simplicity, we  ap¬ 
ply  formula  (13)  to  different  leg  recirculation  schemes 
for  SLIP  without  pitching  dynamics,  i.e.  d  =  0  and 
no  0  dependence.  The  investigation  of  leg  recircula¬ 
tion  strategies  of  SLIP  with  pitching  dynamics  will  be 
published  elsewhere.  For  SLIP  without  pitching,  the 
determinant  at  the  fixed  point  simplifies  to 


det(Ilr2(i2))  =  1-1-  (14) 

\/^  -  ^0  j-dipitTp)  +  Zodi„(l>{iTD))  +  Zq 
-\/l  -  Z^di<j){iTD)  -  Zo 


3.1.3  Analysis  of  Recirculation  Strategies 

Consider  now  the  following  family  of  leg  recirculation 
schemes 

(j){i)  =  k  arccos(zo)  +  a{i  —  Izq)  :  k,l  >  0 
with  a{x)  =  0  for  a;  <  0  (15) 

k  =  I  means  that  recirculation  starts  at  the  leg  liftoff 
angle;  if  /  =  1,  a  certain  angular  trajectory  is  specified 
starting  at  apex  (see  the  leg  retraction  scheme  in  [10]). 
The  application  of  the  partial  return  map  ^2  o  on 
i2o  in  the  threshold  function  (3)  in  order  to  de¬ 
termine  for  which  parameters  of  the  leg  recirculation 
schemes  (15)  volume  is  preserved  is  relegated  to  the 
appendix.  Instead  we  proceed  by  explicitly  computing 
det(iAr2(a;2)).  With  the  angular  trajectory  (15),  the 
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determinant  becomes 


det(-Dr2(i2))  =  1+ 


-Za/1  -  zla{zo{2  -  1))  +  zo{k  -  1) 


\/l  -  z^a{zo{2  -  1))  +  zo 


Hence  for  different  leg  angle  protocols  we  obtain 


(16) 


1.  Constant  leg  touchdown  angle  protocol:  k  =  I  = 
0,  a  =  2tt  —  P  ^  det(Zlr2(x2))  =  0.  The  two- 
dimensional  Jacobian  has  rank  one,  and  the  re¬ 
turn  map  becomes  one-dimensional.  In  [8]  this 
return  map  was  parametrized  by  apex  height, 
whereas  in  [9]  the  angle  of  the  touchdown  velocity 
was  chosen.  No  information  about  the  behavior 
of  this  lower-dimensional  return  map  can  be  ob¬ 
tained  from  this  argument. 

2.  Leg  retraction  [10]:  k  =  0,  I  =  1,  and  a{i  —  Zq)  = 

a  A  +  —  Zq)  where  a  a  is  a  constant  angle  and 

uj  =  is  a  constant  dimensionless  angular 

velocity.  Then  again  det(iJr2(i2))  =  0  and  the 
behavior  of  the  remaining  one-dimensional  return 
map  cannot  be  determined  from  this  argument. 


3.  Leg  recirculation  (starting  at  leg  liftoff):  I  =  0, 
A:  >  0,  and  a{t)  =  aA  +  ojt.  This  exemplifies 
the  fast  rotation  phase  of  the  open  loop  policy 
used  by  RHex  [1],  although  a  full  analysis  is  be¬ 
yond  the  scope  of  the  present  paper.  Then  for 
Zo  >  0  and  zg  <  1  (necessary  conditions  for  a 
pseudosymmetric  flight  phase): 


det(Dr2(x2)) 


zo(l  -  k) 


Zg  Co 

\/l  -  zl 

<  1  : 

Q<k<l 

1  : 

k=l 

>  1  : 

k>  1 

(17) 


In  order  to  illustrate  the  predictive  power  of  eq. 
(17),  we  numerically  approximate  the  determinant 
det{DTZ{x))  of  the  full  return  map  for  fixed  SLIP 
parameters  E  =  ^  =  2.1,  7  =  =  13,  and 

fixed  recirculation  parameters  (xa  =  tt,  o;  =  14  for 
different  k  G  {1/6, 0.5, 1,  2,  3.3}.  Here,  E  is  the  to¬ 
tal  energy  of  the  system  and  the  spring  potential  is 
V{ri)  =  (k/2)(?7  —  ?7o)^.  We  then  compare  these  values 
to  the  values  of  the  determinant  obtained  by  inserting 
the  numerically  determined  fixed  points  X2  =  (zo  Zq)^ 
into  eq.  (17).  The  determinants  obtained  in  those  two 
different  ways  are  plotted  in  Fig.  2a  and  agree  to  a 
high  precision  (|  det(Z17^(a;))  — det(Zlr}(i2))|  <  10“”^). 
Barring  an  improbable  numerical  cancellation  between 
stance  and  flight  phase  dynamics,  this  also  demon¬ 
strates  that  the  SLIP’s  stance  phase  is  volume  pre¬ 


serving.^  In  Figs.  2b-d  iterations  of  the  return  map 
in  (zq  zo)-space  are  shown  for  k  G  (1/6, 1,3.3}  and 
initial  conditions  off  the  fixed  point.  The  eigenvalues 
are  complex  conjugate  pairs  in  all  three  cases.  For 
fc  =  1/6  the  trajectory  spirals  towards  the  fixed  point, 
as  expected  from  a  stable  fixed  point  (Fig.  2c),  for 
fc  =  1  the  trajectory  is  a  deformed  circle  around  the 
fixed  point, ^  indicating  neutral  stability  (Fig.  2d),  and 
for  fc  =  3.3  the  trajectory  spirals  away  from  the  fixed 
point,  indicating  instability  (Fig.  2b). 

3.2  Lateral  leg-spring  model 

The  lateral  leg-spring  (LLS)  was  introduced  in  [6]. 
We  focus  here  on  the  three-degree-of-freedom  version 
with  pairs  of  ‘virtual’  elastic  legs  in  intermittent  con¬ 
tact  with  the  ground.  A  full  stride  consists  of  two 
stance  phases:  a  phase  where  the  first  elastic  leg  piv¬ 
ots  around  a  “foothold”  on  one  side  of  the  rigid  body, 
followed  by  a  phase  where  the  second  elastic  leg  piv¬ 
ots  around  a  “foothold”  on  the  opposite  side.  See  [6] 
for  details.  The  equations  of  motion  of  both  stance 
phases  can  be  cast  into  the  form  (6)  (with  5  =  0). 
Hence  the  stance  phases  from  leg  touchdown  to  liftoff 
are  volume  preserving.  They  are  related  by  a  transi¬ 
tion  mapping  which  maps  the  state  at  liftoff  of  the 
1.  leg  to  the  state  at  touchdown  of  the  2.  leg,  and 
an  analogous  map  Tp.  Thus  the  return  map  reads 
TZ  =  Tp  or\oTp  o  rf . 

The  dynamics  of  the  LLS  model  can  be  described  by 
four  state  variables  (v,  S,  6,  oj),  where  v  is  the  center  of 
mass  speed,  S  is  the  angle  between  the  body  axis  and 
the  mass  center  velocity  vector,  0  is  the  angle  between 
the  body  axis  and  an  inertial  frame  and  oj  =  9.  In 
[6]  these  four  variables  are  augmented  by  two  fixed 
parameters:  /3,  the  leg  touchdown  angle  with  respect 
to  the  body  axis,  and  Ig,  the  relaxed  leg  length;  and 
from  the  values  of  these  six  quantities  at  liftoff  one  can 
find  the  initial  data  for  the  next  stance  phase. 

In  these  variables,  the  transition  mapping  Tp  (omit¬ 
ting  Ig)  reads 

= 

PV  (18) 

where  n  stands  for  the  nth  stance  phase,  LO  for  liftoff, 
and  P  is  held  constant  for  all  stance  phases.  If,  as  im¬ 
plicitly  assumed  for  the  SLIP  treated  above,  P  (and/or 
Ig)  are  regarded  as  state  variables  rather  than  param¬ 
eters,  then  since  they  are  ‘reset’  to  fixed  values  at 

^This  is  not  true  for  approximations  to  the  stance  phase 
dynamics  which  violate  the  time  reversing  symmetry  Si. 

^Despite  the  fact  that  the  determinant  away  from  the  fixed 
point  is  not  1  in  this  example,  volume  seems  to  be  preserved  in 
a  finite  region  around  the  fixed  point. 
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Figure  2:  a)  Comparison  of  the  numerically  computed  determinant  det{DTZ(x))  (+)  of  the  return  map  Jacobian 
to  to  the  determinant  det(Ur2(i2))  (o)  obtained  by  using  the  numerically  determined  fixed  points  in  eq.  (17). 
b)-d)  Trajectories  around  a  fixed  point.  Because  of  the  slow  convergence,  only  every  9th  iteration  in  plot  b)  and 
every  5th  iteration  in  plot  c)  is  shown. 


touchdown,  independent  of  their  values  at  liftoff,  the 
transition  mapping  has  rank  four,  volume  is  not 
preserved,  and  no  deductions  can  be  made  regarding 
the  reduced  four-dimensional  map.  Restoring  a  non¬ 
trivial  dynamical  role  to  the  variable  /?,  for  example, 
via  a  leg  swing  feedback  strategy  similar  to  (15),  could 
lead  to  a  non-degenerate  mapping. 

4  Conclusions 

In  this  paper  we  used  the  example  of  the  SLIP  locomo¬ 
tion  model  to  show  how  factored  analysis  of  the  return 
map  may  be  a  useful  new  tool  in  the  stability  analy¬ 
sis  of  hybrid  Lagrangian  systems.  Specifically,  we  ob¬ 
tained  a  necessary  condition  for  the  asymptotic  stabil¬ 
ity  of  SLIP  in  the  presence  of  a  leg  recirculation  strat¬ 
egy  relevant  to  the  operation  of  the  robot  RHex  [1]. 
This  condition  is  formulated  in  Sec.  3.1.3  for  a  particu¬ 
lar  family  of  leg  recirculation  strategies  as  an  exact  al¬ 


gebraic  expression  despite  the  non-integrability  of  the 
SLIP  system.  Hence  leg  recirculation  strategies  that 
violate  the  above  condition  can  be  discarded  without 
recourse  to  cumbersome  numerical  simulations.  Ap¬ 
plication  of  this  formalism  to  the  robot  RHex  requires 
a  more  elaborate  parametrization  of  leg  recirculation 
schemes  modeled  after  RHex’s  open  loop  controller. 
This  analysis  can  provide  for  the  first  time  a  partial 
explanation  for  the  surprising  self-stable  behavior  ob¬ 
served  empirically  in  RHex.  It  also  paves  the  way  for  a 
more  principled  investigation  of  detailed,  biologically 
motivated  leg  placement  strategies  in  the  LLS  model 
[6]  which  captures  many  aspects  of  cockroach  locomo¬ 
tion  [19]. 
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I  det(Ilr2)|  =  1,  as  was  explicitly  derived  in  (17).  For 
other  values  of  k  and  I  this  does  not  in  general  re¬ 
duce  to  (22),  although  we  have  not  ruled  out  that  for 
specific  values  of  k  and  I  and  a  specific  form  of  a  the 
original  threshold  function  (22)  is  recovered. 


Appendix:  Invariance  of  the 

threshold  equation  for  SLIP 
without  pitching 

In  this  appendix  we  show  that  the  invariance  of 
the  threshold  equation  (11)  under  S2  o  :  X20 
S2  °  /2^°(52o)  with  X20  =  (So  ^0)^  corresponds  to 
|det(Zlr2)|  =  1  for  the  leg  recirculation  family  (15). 
First  we  observe  that 

/l(i2o)  =  (19) 

52  0/I(i2o)  =  (20) 

f\oS2ofl(i20)  =  (21) 
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